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' Abstract 
(N ■ 

J>" I We use renormalization-group methods in effective field theory to improve the theoret- 

ical prediction for the cross section for Higgs-boson production at hadron colliders. In 
^Sj . addition to soft-gluon resummation at N^LL, we also resum enhanced contributions of 

I the form (CavtOs)", which arise in the analytic continuation of the gluon form factor to 

On ' time-like momentum transfer. This resummation is achieved by evaluating the match- 

. ing corrections arising at the Higgs-boson mass scale at a time-like renormalization point 

I //^ < 0, followed by renormalization-group evolution to fi'^ > 0. We match our resummed 

result to NNLO fixed-order perturbation theory and give numerical predictions for the 
total production cross section as a function of the Higgs-boson mass. Resummation ef- 
fects are significant even at NNLO, where our improved predictions for the cross sections 
I at the Tevatron and the LHC exceed the fixed-order predictions by about 13% and 8%, 

respectively, for niH = 120 GeV. We also discuss the application of our technique to 
other time-like processes such as Drell-Yan production, e^e~ hadrons, and hadronic 
decays of the Higgs boson. 
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1 Introduction 



The search for the Higgs boson is of highest priority in the experimental programs at the 
Tevatron and the LHC. A large effort is thus made to obtain precise theoretical predictions 
for the corresponding production cross sections. At hadron colliders the dominant production 
channel is the gluon fusion process through a top-quark loop. The total cross section has been 
calculated in the heavy top-quark limit up to ncxt-to-ncxt-to-leading order (NNLO) in QCD 
[1-6], and fully differential predictions are available at the same order for the decays of the 
Higgs boson into two photons [7] and four leptons [8, 9]. The exact dependence of the total 
cross section on the top-quark and Higgs-boson masses is known at NLO [10, 11]. 

The perturbative corrections to the total cross section turn out to be surprisingly large: for 
a light Higgs boson, the NNLO fC-factor is around 2, and a scale uncertainty of ±(10-15)% 
remains even at this order. The theoretical prediction has been refined using soft-gluon re- 
summation, which has been implemented at NNLL order [12] and recently even at the N^LL 
level [13-17]. The soft-gluon resummation reduces the scale dependence, but the large K- 
factor remains almost unchanged. Indeed, it is not obvious why the cross section should be 
dominated by soft-gluon radiation: given the large center-of-mass energy of the LHC, there is 
plenty of phase space available for hard radiation. 

In a recent paper [18], we have shown that the large K-factor is mostly due to terms of the 
form (Cattcks)" in the perturbative series, which arise in the analytic continuation of (double) 
logarithmic terms in the gluon form factor from space-like to time-like kinematics, InQ^ 
In — in. Being related to Sudakov logarithms, these "vr^-enhanced" contributions can be 
resummed [19-22]. Effective field-theory methods provide a particularly simple framework for 
performing this resummation by implementing matching calculations at time-like momentum 
transfer and extending renormalization-group (RG) evolution into the complex momentum 
plane [18]. At first sight it might appear unsystematic to resum vr^-cnhanccd perturbative 
corrections, which cannot be separated from other numerical coefficients (including tt^ terms 
not associated with analytic continuation) in a parametric way. However, in our RG framework 
this resummation simply corresponds to the proper choice of a particular matching scale and 
as such is unambiguous and physically motivated. The final result for the RG-improved cross 
section follows straightforwardly by applying the usual rules of effective field theory at every 
matching step. Our approach to the resummation of 7r^-enhanced corrections is similar in spirit 
to the analysis of the e+e^ hadrons cross section using a running coupling evaluated at 
time-like momentum transfer discussed a long time ago in [23], and to the "contour- improved 
perturbation theory" introduced in the analysis of hadronic r decays in [24, 25]. It can be 
viewed as an extension of these methods to problems with Sudakov double logarithms. 

In the present paper, we resum both threshold logarithms from soft-gluon emission and the 
TT^-enhanccd terms using the momentum-space formalism developed in [26-28]. Our result is 
based on the factorization of the cross section near threshold into a hard and a soft function. 
The resummation is achieved by solving the RG equations for the different parts. In contrast 
to the standard treatment in Mellin space, this approach yields simple analytic results for the 
resummed hard-scattering kernels in momentum space and is free of Landau-pole ambiguities. 
With a phcnomcnological analysis we investigate to what extent the partonic threshold is 
enhanced due to the fall-off of the parton distribution functions (PDFs). We find that at LHC 
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energies and for the relevant values of the Higgs-boson mass the scale of the soft emission is 
not much lower than rriH, so that no numerically large logarithms arise from soft emissions. 
The main numerical effect of RG improvement is thus due to the resummation of the (C^Tras)"- 
terms in the virtual corrections. In our RG framework, this resummation is accomplished by 
evaluating the hard matching corrections at a scale /i^ — —nijj — ie instead of the conventional 
choice 111 — +mjj. 

We begin our analysis with a brief review of the fixed-order results for the total cross section 
and study to which extent the cross section is dominated by the leading singular terms near the 
partonic threshold. We then discuss the factorization properties of the hard-scattering kernels 
in the threshold region and derive the formulas for the RG resummation of large perturbative 
corrections in momentum space. After determining the default values of the matching scales, 
we present a detailed phenomenological analysis and make predictions for the Higgs-boson 
production cross sections at the Tevatron and the LHC. Compared with previous studies, 
we find significantly faster convergence and improved stability of the perturbative expansion. 
We finally comment on applications of RG-improved perturbation theory to other time-like 
processes, such as Drell-Yan production, the e~e~ — > hadrons cross section, and the total 
hadronic Higgs-boson decay rate. In particular, we explain why the latter two processes do not 
contain 7r^-enhanced corrections of the type present in Drell-Yan or Higgs-boson production. 



2 Fixed-order results 



We consider the production of a Higgs boson in hadron-hadron collisions at center-of-mass 
energy ^/s. The total cross section can be written as 



o- = (7o^ / —Cij{z,mt,mH,lif)fij{r/z,Hf), 
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arc the effective parton luminosities, and are hard-scattering kernels, which are known to 
NNLO in perturbation theory [4-6] . The quantity 
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denotes the Born- level cross section in units of the gluon-gluon luminosity fggir^jif). The 
function A{xq) results from a quark loop connecting two gluons with the Higgs boson. It 
approaches 1 for Xq ^ oo and vanishes proportional to for x^ 0. It follows that Higgs- 
boson production is predominantly mediated by a top-quark loop, while the contributions from 
lighter fermions are strongly suppressed. We include radiative corrections in the heavy top- 
quark limit, i.e., we will only keep logarithmic top- mass dependence in the hard-scattering 
kernels Cij{z,mt,mH, ^f)- For not too heavy Higgs-boson masses the terms suppressed by 
powers of the top-quark mass are numerically very small. Leaving them out greatly simplifies 
the calculation, since one can then use an effective Lagrangian obtained after integrating out 
the top quark. 

Because they are suppressed by Xq, the only numerically relevant correction due to lighter 
fermions is the bottom-quark loop contribution. Its main effect is due to its interference with 
the top-quark loop, which is well approximated by writing | ^^^^(xq)^ ~ (1 — Cb) |74(xt)p with 



While using the pole mass is appropriate for the top quark, the virtual 6-quarks in Higgs-boson 
production are far off their mass shell, and one should thus use the MS quark mass at the 
Higgs mass scale when evaluating the bottom-quark contribution. At uih — 120 GeV we take 
mij ~ 2.8 GeV, and the presence of the bottom-quark loop term in (3) reduces the cross section 
by 6.5%. For comparison, the above approximate treatment would yield = 7.0%. Since the 
6-quark contribution scales like l/m\, it becomes smaller for higher Higgs-boson masses. 

To validate these approximations numerically we have used the computer code [29] , which 
includes the exact quark-mass dependence at NLO [10, 11]. For the range 120 GeV < mn < 
300 GeV we find that the full NLO fixed-order result is about 1% lower than what is obtained 
with the above Born-level treatment of finite top-mass effects. The difference is negligible 
compared to other uncertainties. Also, using the same code we find that the inclusion of the 
6-quark loop decreases the NLO cross section by about 6% at rriH = 120 GeV and about 2% 
at ruH — 300 GeV, in good agreement with the approximate lowest-order treatment described 
above. 

The variable z = m\/s in (1) measures the ratio of the Higgs-boson mass to the parton- 
parton center-of-mass energy v^. The limit 2; — > 1 is referred to as the "partonic threshold 
region". This is the region near Born kinematics, in which the colliding partons have just 
enough energy to produce the Higgs boson. It is an empirical fact that in many cases this 
region gives the dominant contributions to the cross section. In Section 3 we will resum these 
contributions to all orders in perturbation theory. 

The sum in (1) extends over all possible combinations of initial partons, but only Cgg 
contributes at the Born level and contains the leading singular terms in the limit 2; — > 1. We 
split off these terms by writing 



where the second piece does not contain singular distributions for 2; — > 1. The explicit expres- 
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Cgg{z, mt, mH,IJ'f) = C{z, rut, mH,IJ'f) + C^fiz, mt, mn, A*/) , 
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sion for the leading singular terms through 0{a'^) is [2-6] 
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where we have defined the distributions 

Dn{z) = 
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The reason for including a factor I/2; in the argument of the logarithm was explained in [28] . 
The remaining contributions to the hard-scattering kernels are free of singular distributions. 
At NLO they read 
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In these expressions = Q;s(/ij), and wc do not distinguish between the factorization and 
renormalization scales. Note that the coefficients of the ln[m^(l — z)^/ii'jz] terms are propor- 
tional to the Altarelli-Parisi sphtting functions. 

To visualize the numerical importance of the leading singular terms for Higgs-boson pro- 
duction, we compare the contributions from these terms with the complete fixed-order results 
in Figure 1. Throughout our analysis we use MSTW2008NNLO PDFs [30] and the associ- 
ated normalization Q;s(m|) = 0.1171 ± 0.0036 of the running coupling constant, unless noted 
otherwise. We use three-loop running and Uf — 5 fight quark flavors. The figure shows that 
the complete fixed-order results are well approximated by the leading singular terms. Taking 
"^H = l^f = 120 GeV as an example, the leading singular terms contribute 96% (94%) of the 
NLO (NNLO) cross section at the Tcvatron, and 90% (86%) of the NLO (NNLO) cross section 
at the LHC. More specifically for the LHC this means that 82% (74%) of the NLO (NNLO) 
correction term are captured by the coefficient C in (7). Note also that only —1% (—8%) 
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Figure 1: Comparison of the complete fixed-order results (solid lines) and the contributions 
from the leading singular terms (dashed lines) to the total cross sections for Higgs-boson 
production at the Tevatron (left) and the LHC (right). We set /x/ = mn- Darker lines 
represent higher orders in as- 



of the NLO (NNLO) correction to the cross section are due to parton production channels 
different from gg H. 

In [28] we have investigated for the case of Drell-Yan production the question to what 
extent the dominance of the leading singular terms can be justified based on the strong fall-off 
of the parton luminosities. In the present case, setting /i/ = 120 GeV for example, we find 
that ffgg{y,fJ'f) oc y~"- with a f» 2.5 for y < 0.05, and ffgg{y,fJ'f) oc (1 — y^ with b ^ 14.5 for 
y > 0.3. Due to this strong fall-off, the integral in (1) is dominated by z values near r. For r 
values exceeding 0.3, the partonic threshold contributions would be enhanced by logarithms of 
b ~ 14.5. However, even at the Tevatron the center-of-mass energy is so high that r < 0.02 for 
Higgs-boson masses below 300 GeV. In this region the cross section (1) is well approximated 
by the simple formula [28] 

~ C^Born / dz z"''^ C {z , TTlt, TTIh , fi f) ] 0"Bom = CTq #99 (^, /W/) , (10) 

Jo 

with a — 1 ~ 1.5. Since the weight function z"""^ is not strongly peaked near z = 1, the 
threshold dominance cannot be explained parametrically in this case. Indeed, we will see later 
that threshold resummation alone has a very minor effect on the predictions for the cross 
section. As a side remark, we note that (10) implies the scaling a oc Tn^^^"' !^ i^h^- 

Let us now discuss in more detail the different momentum regions that contribute to the 
Higgs-boson production cross section. For a not too heavy Higgs boson, the gluon-gluon fusion 
process gg —>■ H is well approximated by the effective local interaction [31-35] 

£eff = aiml /i^) - G,.,, Gr , (11) 

V lZtc 

where v ~ 246 GeV is the Higgs vacuum expectation value, and /i denotes the scale at which 
the local two-gluon operator is renormalized. The short-distance coefficient Ct is known up to 
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NNNLO [36, 37]. To NNLO, the result reads [38, 39] 
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The production cross section is related to the discontinuity of the product of two such effective 
vertices. As explained earlier, it is a good approximation to keep the exact dependence on 
the top-quark mass in the Born- level cross section ctq in (1), but to employ the effective local 
interaction (11) for the analysis of higher-order perturbative corrections. 

Once the top quark has been integrated out, the hard-scattering kernels receive contribu- 
tions associated with two different scales: a "hard" scale fi\ ~ set by the mass of the 
Higgs boson, and a "soft" scale /i^ ~ s(l — 2;)^ = m|^(l — z)'^/z, where s(l — 2^/2 = (j9^)max 
is determined by the maximum available transverse momentum. The presence of these two 
scales is apparent from the structure of the logarithms in the fixed-order terms in (7) and (9), 
but it can also be derived more rigorously using the method of regions [40] . The short-distance 
coefficient C in (7) can be factorized as 

C{z,mt,mH,l^f) = [Ct{m^,i^j)Y H{mH,l^}) S{s{l - zf,i^j) , (13) 

and in this way the scale separation becomes explicit. The derivation of this factorization 
theorem using effective field theory proceeds in analogy with the discussion for the Drell-Yan 
case presented in [28] and will not be repeated here in detail. The formula results from a 
sequence of matching steps illustrated in Figure 2. The Standard Model with six quark flavors 
is first matched onto a five-fiavor theory by integrating out the heavy top quark. The Wilson 
coefficient arising in this step is Ct- In the next step, the five-flavor Standard Model is matched 
onto soft-collinear effective theory (SCET) [41-44] containing soft degrees of freedom along 
with two types of hard-coUinear fields aligned with the directions of the particle beams. The 
corresponding Wilson coefficient is H. In the third step, this theory is matched onto another 
version of SCET, in which the soft modes are integrated out and the hard-coUinear modes are 
replaced by collincar fields of lower virtuality. The soft function 5* is the matching coefficient 
arising in this step. The remaining low-energy matrix element is then identified with the 
parton luminosity function ffgg defined in (2). The calculation of the components Ct, H, and 
S at any order in perturbation theory is much simpler than the calculation of the Higgs-boson 
production cross section at the same order. The factorization formula (13) thus provides an 
approximation to the cross section that requires a minimal amount of calculational work. The 
all-order resummation of the partonic threshold logarithms ("soft-gluon resummation" ) and 
of other, "7r^-enhanced" terms [18] is then achieved by solving RG equations. 

SCET provides field-theoretic definitions of the factors H and S in the factorization for- 
mula. The hard function H is the square of the on-shell gluon form factor evaluated at time-like 
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Figure 2: Sequence of matching steps and associated effective theories leading to the factor- 
ization theorem (13). 



SM 




n/ = 6 






SCET 

he, he, s 



momentum transfer = m^, and with infrared divergences subtracted using the MS scheme 
[17, 26, 28]: 



H{mH,fi ) = ICsi-rriH - ie,fi 



(14) 



On a technical level, the function Cs appears as a Wilson coefficient in the matching of the 
two-gluon operator in (11) onto an operator in SCET, in which all hard modes have been 
integrated out. This matching takes the form 



(15) 



where = —q^ is (minus) the square of the total momentum carried by the operator. The 
fields and ^1^'" are effective, gauge-invariant gluon fields in SCET [45]. They describe 
gluons propagating along the two light-like directions n, n defined by the colliding hadrons. 
The two-loop expression for the Wilson coefficient Cs can be extracted from the results of 
[46]. We write its perturbative series in the form 
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where L = \n[{—m'jj — ie)//i^]. The one- and two-loop coefficients read 
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The soft function S in (13) is defined in terms of the Fourier transform of a vacuum 
expectation value of a Wilson loop in the adjoint representation of SU{Nc). In SCET is 
arises after the decoupling of soft gluons from the hard-coUinear and anti-hard-collinear fields 
describing the partons originating from the colliding beam particles [28]. The soft function 
in the case of Higgs-boson production is closely related to an analogous function entering 
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the Drell-Yan cross section [17, 28]. At two- loop order (but not beyond) the two quantities 
coincide after a simple replacement of color factors. In the notation of the second reference, 
we have 

S{s{l - z)\ ix)) = VI M/Higgs(s(l - zf, ix)) 

(18) 

The explicit form of the result can be derived using formulas compiled in Appendix B of [28] . 

When one inserts the two-loop expressions for the various component functions into (13) 
and expands the product to 0{a1), one recovers the expression given in (7). In the following 
section we will discuss how improved perturbative expressions for the component functions 
can be obtained by solving RG evolution equations with appropriate boundary conditions. In 
this way one avoids perturbative logarithms arising when the factorization scale ///is chosen 
different from the characteristic scales m^, m^/, or -\/I(l — z). Even though these logarithms 
are not particularly large, their resummation has the effect of improving the stability of the 
prediction with respect to scale variations. More importantly, however, we will also be able to 
resum the vr^-cnhanccd terms in the perturbative expansion related to the time- like kinematics 
of the Higgs-boson production process. They have been shown to be responsible for the bulk 
of the large X-factors arising in calculations of the Higgs-production cross sections at the 
Tevatron and the LHC [18]. 



3 Renormalizat ion- group analysis and resummation 

Our formalism for the resummation of large perturbative corrections in Higgs-boson production 
is based on effective field theory and follows closely our previous analyses of DIS at large x 
[26, 27] and Drell-Yan production [28]. The two key steps of the approach are deriving a 
factorization formula such as (13) valid near the partonic threshold z ^ 1, and then using the 
RG directly in momentum space to resum logarithms arising from ratios of the different scales. 
We stress that the final, RG-improved formula for the cross section follows unambiguously by 
applying the rules of effective field theory at each step of the derivation. 

The Wilson coefficient Ct appearing when the top quark is integrated out satisfies the RG 
equation 

^ a(m^/.^) = f(«.)G*(m^/.^), with l\cXs)^ocl4r^- (19) 



7 -i — i-v t't" / I \ — B/ ~v\ 1 ' r- J 1 / \ — a/ — s J •) 

a In jj, das (^i 

The fact that the anomalous dimension is related to the QCD /3-function [34, 47] is not 
surprising, since the two-gluon operator in (11) is proportional to the Yang- Mills Lagrangian. 
The evolution equation can be integrated in closed form and leads to 

(5[as{ii.i))/al{ii,i) 

where ~ rrit is the matching scale at which the top quark is integrated out. 
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The Wilson coefficient Cg arising when hard, virtual quantum corrections to the effective 

two-gluon vertex (11) are integrated out obeys an evolution equation reflecting the renor- 
malization properties of the effective two-gluon SCET operator on the right-hand side of the 
matching relation (15). It reads 
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Cs{-ml-ie,f,'), (21) 



where T-^^^^ is the cusp anomalous dimension of Wilson lines with light-like segments in the 
adjoint representation of SU (Nc). It controls the leading Sudakov double logarithms contained 
in Cs and is known to three-loop order [48]. The single-logarithmic evolution is controlled 
by the anomalous dimension 7"^, which can be extracted from the infrared divergences of the 
on-shell gluon form factor [26]. Using results from [49] it can be derived to three-loop order. 
We collect the relevant expressions for the expansion coefficients of the anomalous dimensions 
in Appendix A. The general solution to (21) is [51] 



Cs{-mH-ie,n'j)=exp 
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where \Xh is the hard matching scale. We have introduced the definitions 
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and similarly for the function a^s. The perturbative expansions of these functions obtained 
at NNLO in RG- improved perturbation theory can be found in the Appendix of [27]. 

The naive choice /i| ~ m'jj of the hard matching scale gives rise to large vr^ terms in the 
matching condition (16), which arise since = ln^[(— — ie)//x^] ~ — tt^ and render the 
perturbative expansion of the hard function H in (14) unstable. We have shown in [18] that 
these TT^-enhanced terms are to a large extent responsible for the poor perturbative behavior of 
fixed-order predictions for the Higgs-boson production cross sections at hadron colliders. We 
can exploit the fact that the solution (22) is formally independent of the hard matching scale 
to avoid the large tt^ terms in the matching condition by a proper choice of the matching scale. 
To this end we set /x^ ~ ~ so that ln[(— — ie)//x^] remains a small parameter. 

The TT^-enhanced terms are then resummed to all orders in perturbation theory and appear 
in the functions S and or in the exponent in (22). With this choice, relation (22) involves the 
running coupling Q;s(/i^) evaluated at negative argument. The definition P{as) — das / din /i 
of the QCD /3-function implies that 

= (24) 
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where here and below as{—ii'^) is to be understood with a —ie prescription. This relation 
allows us to define the running coupling at time-like argument in terms of that at space-like 
momentum transfer. At NLO we obtain 

where we count a(/i^) = /3oQ;s(Ai^)/4 as an 0{1) parameter. It is important that (24) is 
independent of the path. For example, the evolution of the coupling can be performed on a 
circle with fixed radius in the complex momentum plane, thereby avoiding the region near the 
origin, where perturbation theory breaks down. Note that the perturbation-theory coupling 
q;s(//^) is analytic in the complex /i^-plane with a cut along the negative real axis. It has 
an unphysical Landau pole at /x^ = A|^, which is of no concern to our discussion since we 
are interested in very large values. In practice, we obtain as{fJ^1) for < by simply 
evaluating the three- loop running coupling at negative values of its argument. 

The soft Wilson loop Wniggs in (18) obeys an integro-differential evolution equation, which 
is analogous to that for the soft function in Drell-Yan production discussed in [28] . The general 
solution to this equation can be obtained using a Laplace transformation [26] . It can be written 
with the help of an associated function Sniggs, which is given by the Laplace transform of the 
soft Wilson loop at a matching scale /is- The solution is then obtained from 



^ W^Higgs(c^', = exp [-^S{iil ii)) + 2a^w{iil n))] SHiggs(9^, lA) [j^j , (26) 

where (9^ denotes a derivative with respect to an auxiliary parameter r/, which is then set to 
77 = 2ar{fJ.% fJ.'j-) ■ As written above, the solution is valid as long as 77 > 0. From the RG 
invariance of the Higgs-boson production cross section one can derive a relation between the 
anomalous dimension entering in the above solution and the anomalous dimensions of the 
remaining components in the factorization formula for the cross section [28] . It reads 

7^ = ^ + 7* + 7^ + 27^, (27) 
as 

where 27^ is coefficient of the 5{1 — x) term in the Altarelli-Parisi splitting function Pg^g{x). 
The three-loop expression for this quantity was obtained in [48] , and we collect the correspond- 
ing expansion coefficients in Appendix A. At two-loop order, relation (18) implies that the 
associated soft function sniggs is obtained from that in the Drell-Yan case by the replacement 
Cp —>■ Ca- This gives 

^mggs(i^, /x^) = 1 + ^ Ca (21^ + y ) + ' {Cl Wa + CaT^uj W,) , (28) 
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The result (28) agrees with a corresponding expressions derived in [17]. 

Putting everything together, we arrive at our final formula for the RG-improved expression 
for the hard-scattering coefficient in (7). It can be written in the form 



C{z,mt,mH,l^f) ^ [Ct(mf,i^t)Y \Csi-mH - ie, U{mH, Ht, l^h, l^s, l^f) 

-2'yEV 



(30) 



where 



U{mH,l^ul^h,l^s,l^f) 



q;2(^2; 



-m 



H 



It 



r(2r^) 



(31) 



X 



Apart from the factor containing the /5-function, which is related to the evolution of the 
two-gluon operator in (11), and the ratio of running couplings, which compensates the scale 
dependence of the Born-level cross section ctq in (1), this result is of the same form as the 
corresponding expression arising in Drell-Yan production and given in equations (50) and 
(51) of [28]. Some comments on the effect of the resummation of 7r^-enhanced terms in the 
Drell-Yan case will be made in Section 6.1. 

It is instructive to consider the special limit in which all matching scales are set equal to a 
common scale /^/ ~ ttih, while = —^'j — ie is still chosen in the time-like region. We then 
obtain [18] 



]nU{mH,fJ^f,-ifJ^f,fJ^f,fJ^f) 



^ 



2131 W{m]j) 



+ 



rf 



[2a arctan(a) — ln(l -|- a^)] 
ln(l + a^) 







(32) 



+ ^ [4arctan2(a) - ln2(l + a^)] + 0(a,)| , 



where a = a{m%). Note that the result is ///-independent at this order. The expression for 



the evolution function simplifies considerably if we treat a(m|^) 
Us- Using the fact that 7q = 0, we then find 



0.2 as a parameter of order 



lnU{mH,l^f,-il^f,l^f,l^f) 



(33) 



This result makes explicit that the Tr^-enhanced corrections are terms of the form (CavtcIs)" 
in perturbation theory and exponentiate at leading order. Numerically, setting Hf — ttih — 
120 GeV we obtain InC/ = {0.558,0.560,0.561} at LO, NLO, and NNLO from the exact 
expression for the evolution function derived from (31), indicating that the leading-order terms 
give by far the dominant effect after RG improvement. The analytical expressions (32) and (33) 



11 



provide accurate approximations to the exact results. The first equation gives InU — 0.557, 
while the second one yields \nU = 0.565. The close agreement of these two numbers shows 
that the running of the coupling between fij- and —fij- is a minor effect compared with the 
evolution driven by the anomalous dimension of the effective two-gluon operator in (15). 

The RG-improved prediction for the leading singular contributions to the Higgs-boson 
production cross section is obtained by integrating the above expression for the resummed 
kernel C with the gluon-gluon luminosity function ffyy, see (1). In order to also account for the 
remaining contributions to the cross section, we add to this result the fixed-order contributions 
arising from the non-singular terms in the hard-scattering kernels, which at NLO have been 
compiled in (9). In the momentum-space approach the subtractions required to avoid double 
counting of the resummed terms are straightforwardly implemented as [27, 28] 

(34) 

We have used this prescription to calculate the fixed-order expressions for the two terms on 
the right-hand side of (6). Note that only after this matching step the cross section is formally 
independent of the factorization scale iif. 

Traditionally, threshold resummation is performed in moment space rather than momen- 
tum space. For the case at hand, one considers moments of the cross section in r = nijj/s at 
fixed niH- 

[ drr^-^a. (35) 



„RGI ^resummed 


1 / ^fixed order 


^resummed 




a —a 


+ [a 








|J■t,^J.h,^J■s,^J.f \ 







After this MeUin transformation, the convolution integrals in (1) and (2) reduce to products 
of moments. For the gluon contribution to the cross section, one has 

(Tjv = (To Cjv+1 (mt , ruH, fif) fj+i M Jn+I M ■ (36) 

Having the analytical result (26) for the RG equation of the soft function at hand allows us 
to work directly in momentum space. However, to compare to the results obtained using 
traditional methods, it is instructive to transform our result (30) for the hard-scattering co- 
efficient to moment space. This was discussed in detail in [27] for DIS and [28] for Drell-Yan 
production. The discussion for Higgs production is completely analogous to the Drell-Yan 
case, but for the numerical discussion below, it will be useful to have explicit formulae also 
for the present case. To obtain the moment-space result we note that [50] 

J^z z^-' S{mUl - z)\ f/) = ?Higgs ( In /.^) + ^ (^) > (37) 

with N = e^'^N . Solving the associated RG equation and combining it with the hard function 
we obtain 

CN{mt,mH,fif) = [Ct{m^,i^^)Y \Cs{-mH - ie, 

f rri^ \ -_ f l\ (38) 
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The result has exactly the same structure as (30). Evaluating the derivatives with respect to 

2 2 

Tj produces logarithms In which can be eliminated choosing ii^ — For fixed /is, the 
result (38) has a very simple A'"-dependence, and we can Mellin-invert it analytically using 



1 r 



c+ioo p-'ilEV 



r{2rj) 



The main difference to (30) is an overall factor of y/z, which amounts to a first-order power 
correction in the threshold region 2; — > 1 . The numerical impact of this factor will be discussed 
below. 

The result (38) can be compared to the expressions used in the traditional formulation 
of resummation. For DIS and Drell-Yan production, it has been shown in [27, 28] that the 
two methods give identical results for the threshold-enhanced terms when expanded to any 
fixed order in a^. In these papers, an exact relation between the radiation functions appearing 
in the traditional framework and the anomalous dimensions and Wilson coefficients in the 
effective theory was derived. For the Higgs case, the corresponding relation reads 

g2,.v ^ 2V) = ^w{as) + V lnS-Higgs(0, - -^-^ T^^^ia^) , (40) 

where = Q!s(At), and V = d/dlnii'^ = \/3{as)/2] d/dag. Using this result, wc reproduce the 
perturbative expression for the radiation function D{as) given in [13, 14] up to third order 
in as- 



4 Choice of the matching and factorization scales 

The RG-improved cross section in (34) is formally independent of each of the three matching 
scales /it, fifi, and fig, as well as of the factorization scale /i/ at which the parton densities arc 
evaluated. However, in practice a residual scale dependence remains due to the truncation of 
perturbation theory. It is a standard procedure to take this residual dependence on the scales 
as an estimate of yet unknown higher-order effects. In our analysis we will independently 
vary the various matching scales about their default values, whose determination we will now 
discuss. In the spirit of effective field theory, the matching scales should be chosen such that 
the matching conditions (i.e., the Wilson coefficients evaluated at the matching scales) in (30) 
have well-behaved perturbative expansions. All large corrections are then resummed into the 
evolution function U. 

The characteristic scale of the top-quark loop integrated out in the construction of the 
effective local interaction (11) is the top-quark mass, and we thus take fit = rrit as our default 
choice for the first matching scale. With this choice, the perturbation series for the matching 
coefficient Cf is well behaved. Setting n/ = 5 for the number of light quark flavors, we flnd 

Ct{ml,ml) = 1 + Q.S7bas{ml) + Q.Q22,al{ml) + . . . . (41) 
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Figure 3: Left: Relative contributions to the total cross section arising from the one-loop 
corrections to the soft function function of the soft matching scale /i^, obtained 

with /i/ = rriH = 120 GeV. The numbers on the curves show the corresponding values of 
r. Right: Results for the soft matching scale Hs for different values of r and four different 
Higgs-mass values (see text). The upper set of curves corresponds to convergence criterion I, 
the lower one to criterion II. 



As mentioned earlier, the most naive choice for the hard matching scale, = f^H, does 
not lead to a well-behaved expansion for the hard matching coefficient. We find 

Csi-mjj - ie, m]j) = 1 + 2.749 a,(m^) + (4.844 + 2.071^) a^(m^) + . . . . (42) 

The origin of the large expansion coefficients can be traced back to the fact that the Sudakov 
(double) logarithms contained in the coefficients c„(L) in (16) give rise to vr^ terms when 
we analytically continue L ln(m|^//i^) — ivr. The same happens for the coefficient Cy in 
Drell-Yan production [19, 20] and for other time-like processes [21]. A vastly better behavior 
is obtained when the matching scale is chosen in the time-like region [18]. This gives (all 
arguments are defined with a —ie prescription) 

Cs{-mjj, -m]j) = 1 + 0.393 a,(-m^) - 0.152 a^(-m^) + . . . . (43) 

Note that the values of the strong coupling in the space-like and time-like regions are not very 
different from each other. For instance, setting mn = 120 GeV we find as{—m'jj)/as{m'jj) = 
0.951 -|- 0.213i. It follows that the stark difference between (42) and (43) is not due to the 
evolution of the running coupling between space-like and time-like values of its argument, but 
rather due to the evolution of the effective two-gluon operator (15) driven by its anomalous 
dimension. In our phenomenological analysis we will thus use /i| = —mlj as our default 
choice. Then the vr^-enhanced corrections are resummed into the evolution function U in (30). 
In order to illustrate the significance of this resummation, we will sometimes use the naive 
choice /i^ = m'jj for comparison. 

Let us now discuss the choice of the soft matching scale fig, which is non-trivial since the 
soft function S in (18) depends on the convolution variable z. For the determination of the 
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soft scale we follow the method proposed in [28], i.e., we choose the value of /ig so that the 
perturbative expansion of the soft function exhibits a good convergence after the integration 
over z has been performed. The result thus depends on the process (in particular, on the 
value of the Higgs-boson mass) and on the shape of the gluon distribution function. The 
left panel in Figure 3 shows the relative contributions to the cross section (normalized to 1) 
arising from the one-loop terms in the soft function SRiggs ^ ^ function of /ig. We choose 
fif = rriH = 120 GeV and consider different values of r = mjj/s between 0.00005 and 0.03, 
which is the relevant region for our study. The two scale-setting criteria proposed in [28] are: 

I. Starting from a high scale, determine the value of /ig at which the one-loop correction 
drops below 15%. 

II. Choose the value of Hs for which the one-loop correction is minimized. 

With either choice the two-loop corrections at the corresponding /ig values are very small, 
indicating a good convergence of the perturbative series. The same analysis is repeated for 
different values of the Higgs-boson mass. The resulting values for the soft scale //^ are shown 
in the right panel in Figure 3 for niH = 120, 160, 200, and 240 GeV. Note that the ratio 
fig/rriH is to a good approximation independent of rriH- An analogous scaling behavior was 
also observed in the Drell-Yan case [28]. Below we will vary the soft scale between the two 
choices labeled by nl and f/g in the figure, taking the average of the two prescriptions as our 
default value. In practice the values of r relevant to Higgs-boson production at the Tevatron or 
LHC are so small that the ratio [Is/itih can be considered an 0{1) parameter. This confirms 
our earlier argument stating that there is no parametric reason to perform soft-gluon (or 
threshold) resummation in this case. 

Following common practice, we take A*/ = mn as our default value for the factorization 
scale. From the perspective of effective field theory it would be more natural to choose /// 
at or below the soft matching scale fig, because the parton densities describe the low-energy 
hadronic matrix elements after the hard and soft modes have been integrated out. 

In our phenomenological analysis we vary the different scales independently about their 
default values and add up the corresponding variations of the cross section in quadrature. 
Specifically, we take mt/2 < jit < 2mt, {niH/^y < — /i| < (2mij)^, and < /ig < 
for the matching scales, and mu/^ < fif < 2m// for the factorization scale. The effect of 
varying is so weak that we do not depict it in a plot. The and fig dependences of the 
cross sections obtained at different orders in RG-improved perturbation theory are shown in 
the upper panels of Figure 4. In both cases the scale dependence strongly decreases in higher 
orders and is essentially negligible already at NLO. In the lower panels of Figure 4, we compare 
the iJ,f dependence of the combined cross section to that of the fixed-order cross section. Note 
that the fif dependence of the fixed-order cross section is not significantly improved when 
going from LO to NLO, and that a sizable fif dependence remains even at NNLO. On the 
other hand, the nj dependence after resummation is already small at LO, while the NLO and 
NNLO resummed cross sections only depend very weakly on the factorization scale. 
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Figure 4: Dependence of the resummed cross section for Higgs-boson production at the LHC on 
the scales ^ih (upper left) and yU^ (upper right), and on the factorization scale /i/ (lower left), for 
niH = 120 GeV. The darker curves correspond to higher orders in RG-improved perturbation 
theory. For comparison, we also show the dependence on /xj in fixed-order perturbation theory 
(lower right). The corresponding curves for the Tevatron would look very similar except for 
the overall scale. 



5 Predictions for the cross section 

We now present numerical results for the Higgs-boson production cross sections at the Teva- 
tron and the LHC. To estimate the theoretical uncertainties we combine the various scale 
dependences as described in the previous section. The effect of the uncertainties in the PDFs 
is estimated by scanning over the 30 different sets provided by [30]. The uncertainty in the 
value of the running coupling ^^('^l) = 0.1171 ± 0.0036 introduces an additional error in the 
cross-section predictions of about ±6%. We compare our RG-improved results for the cross 
sections with those obtained in fixed-order perturbation theory. In the latter case we vary the 
factorization and renormalization scales together in the range mj7/2 < fif < 2m h- 

In Figure 5 we show the scale dependence of our predictions for the cross sections at 
different orders in perturbation theory. Note that we use the same PDFs (MSTW2008NNLO) 
in all cases, i.e., we do not switch to LO or NLO distribution functions in the lower-order 
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Figure 5: The fixed-order (left) and RG-improved (rigfit) cross-section predictions including 
perturbative uncertainty bands due to scale variations for tfie Tevatron (upper) and LHC 
(lower plots). Darker bands correspond to higher orders in perturbation theory. 



results. This makes it easier to judge the actual size of the perturbative corrections to the 
hard-scattering kernels. The results obtained after RG improvement show significantly faster 
convergence and reduced scale dependence in higher orders. The NNLO resummed predictions 
have a perturbative uncertainty of less than 3% for both the Tevatron and the LHC, while the 
scale dependence of the NNLO fixed-order results is approximately ±15% for the Tevatron 
and ±10% for the LHC. Numerical values for the cross section at NNLO are shown in Table 1. 
The first error accounts for scale variations, while the second one reflects the uncertainty in the 
PDFs. The additional uncertainty of ±6% due to the value of as{rn\) is not shown explicitly. 
We emphasize that the effect of RG improvement is significant even at NNLO, where the 
resummed cross sections at the Tevatron and the LHC exceed the fixed-order predictions by 
about 13% and 8%, respectively (for rriH = 120 GeV). These differences are as important 
numerically as the differences between the NLO and NNLO resummed results. 

In Figure 6, we show for comparison the results obtained when the PDFs are switched 
according to the order of the calculation. When this is done, the higher-order bands obtained 
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Figure 6: The fixed-order (left) and RG-improved (rigfit) cross-section predictions including 
perturbative uncertainty bands due to scale variations for the Tevatron (upper) and LHC 
(lower plots). In contrast to Figure 5, different PDF sets are used according to the order of 
the calculation. 



after RG improvement are fully contained in the lower-order ones and the i^-factor is close 
to 1, in particular for the LHC.^ In fixed-order calculations it is customary to use PDFs ex- 
tracted from a fit using predictions of the same order. Doing so absorbs universal higher-order 
corrections into the PDFs. Since resummed calculations contain contributions of arbitrarily 
high orders, the optimal PDF choice is less clear. If the same large higher-order corrections 
affect both the observable one tries to predict and the cross sections used to extract the PDFs, 
it would be quite problematic to perform a resummation in one case and not the other. For 
our case, the relevant input quantity is the gluon PDF at low x, which is mostly determined 
by measurements of scaling violations in the DIS structure function, dF2{x,Q'^)/dQ'^ . The 
higher-order corrections associated with the analytic continuation of the time-like gluon form 
factor, which we resum, do not affect the DIS cross section, and so are not universal and 

^For MRST2004 PDFs [52], the iiT-factors after resummation are somewhat larger, K w 1.3 for the LHC, 
see [18]. 
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Table 1: Cross sections (in pb) for different Higgs masses at the Tevatron and tlie LHC, using 
MSTW2008NNLO PDFs. The first error accounts for scale variations, while the second one 
reflects the uncertainty in the PDFs. 



ruH [GeV] 



LHC 

fixed-order NNLO resummed NNLO 



Tevatron 

fixed-order NNLO resummed NNLO 
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cannot simply be absorbed into the PDFs. However, such corrections will be present for the 
Drell-Yan cross section, but smaller by a factor Cf/Ca — 4/9. In Section 6.1, we discuss their 
resummation for this process. Similar effects also appear in jet production processes, and it 
would be interesting to extend our method also to this case. Since jet production involves an 
interplay of time-like and space-like dynamics, the identification of the enhanced contributions 
will be more involved in this case. 

While our numerical values include an uncertainty due to PDFs, it is important to note that 
these uncertainties arc estimated within the theoretical framework adopted when performing 
the PDF extraction. There are examples where newer PDF determinations have led to shifts 
which are larger than the quoted uncertainties. In particular, with the new MSTW2008NNLO 
PDFs the Higgs production cross section at the Tevatron is reduced by 10-15% compared to 
MRST2006NNLO [53], the previous-generation PDF set by the same collaboration. Also, the 
MRST2006NNLO PDFs had lead to an increase of the production cross section at the LHC 
by 10% compared to the result obtained using MRST2004NNLO [52]. For comparison, we 
provide in Table 2 results obtained using CTEQ6.6 PDFs [54]. They are in good agreement 
with the results given in Table 1. Note, however, that the CTEQ6.6 PDFs are obtained from 
a fit to data using NLO cross sections only. Note also that the MSTW2008NNLO gluon PDF 
differs significantly from the one obtained by Alekhin et al. [55, 56]. 

In our predictions we resum logarithmic terms near the partonic threshold as well as 
the TT^-enhanced terms contained in the hard matching coefficient H in (13). It is simple 
to disentangle the two effects: choosing /x^ = rn?jj instead of ii\ = —rrijj switches off the 
resummation of the tt^ terms. With this choice our results are equivalent to what is obtained 
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Table 2: Cross sections (in pb) for different Higgs masses at the Tevatron and tlie LHC, using 
CTEQ6.6 PDFs. Tfie first error accounts for scale variations, while the second one reflects 
the uncertainty in the PDFs. 
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3_o.4 


1 c r+0.4+0.4 
-■-^■^-0.1-0.4 


A nSS+0-012+0.016 
U.UOO_Q 012-0.014 


A 1 Aq+0.002+0.018 

U.1UO_0.000-0.015 
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1 9 1+0.3 
^^■"-1.1-0.3 


1 7+0.3+0.3 
-■-"J- '-0.1-0.4 


U.UOO_o.o09-0.011 


A A70+O.OO2+O.OI5 
U.U/O_0.000-0.012 


280 


11 4+1.0+0.3 
^J^-^-1.0-0.3 


1 9 4+0.3+0.3 
-■-^■^-0.1-0.4 


U.U01_o.007-0.010 


u.uou_o.ooo-o.oii 


300 


1 n c:+0.9+0.3 

i>J.O_l Q_Q 3 


11 4+0.2+0.3 
^^■^-0.1-0.3 


n nzLn+0 006+0.010 

u.u^u_o 006-0.008 


u.u'±o_o 000-0.009 



in standard soft-gluon resummation, albeit performed in momentum space instead of Mellin 
moment space. As seen in Table 3, the main effect of RG improvement is the resummation of 
the TT^-enhanced terms contained in the hard matching coefficient H in (13). The predictions 
for the resummed cross section obtained without resummation of the vr^-enhanced terms are 
quite close to the fixed-order results. This shows once again that soft-gluon resummation is 
not an important effect in the case of Higgs-boson production at Tevatron or LHC energies. 
It confirms our theoretical argument given in conjunction with relation (10) and is also in line 
with our finding that the optimal value of the soft scale is not much lower than hard scale 
set by the Higgs mass. 

Compared to the numerical results obtained in [12] using the traditional moment-space 
formalism, we find smaller threshold resummation effects. Part of the difference is simply 
due to the fact that we use the values obtained with our default scale choices as our default 
values. As a consequence our uncertainty bands are quite asymmetric. In contrast, [12] 
uses the central values of the bands as the default values. However, even after taking this 
trivial difference into account, we observe that the resummation effects found in [12] are larger 
than those we find. While compatible within the assigned uncertainties, we find that the 
resummation effects obtained using the moment formalism of [12] are about twice as large as 
those we find in momentum space. Specifically, evaluating the NNLL moment-space result of 
[12] with MSTW2008NNLO PDFs, and matching to the NNLO fixed-order result, we obtain 
a = 50.4 pb at the LHC for mn = 120 GeV, compared to our result a = 48.5 pb and the 
fixed-order value of 47.6 pb. 

There are three differences between our calculation and what was done in [12]: (i) we go one 
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Table 3: Cross sections (in pb) for mn = 120 GeV. We compare fixed-order results (first 
column) with RG-improved results (remaining three columns) corresponding to standard soft- 
gluon resummation with /x^ = +m%, resummation of 7r^-enhanced terms (//^ = ~"^lf) oiily? 
and the combination of both. The uncertainties are due to scale and PDF variation. 







fixed order 


threshold 


TT^-enhanced 


threshold + tt^ 


LHC 


LO 
NLO 
NNLO 


1 ^ c,+2.4+0.4 
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AK 0+3.0+1.1 
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LO 
NLO 
NNLO 



n 90-1 +0.105+0.018 
U.ZOi_o. 071-0.019 



0.389 



0.650 



hO. 172+0.041 
-0.131-0.044 

n Oni +0.126+0.056 
U.yui_0.i24_o.060 



-0.062+0.023 
-0.046-0.024 
n 7«A+0-0''7+0.045 
U. ( U'±_Q Q2g_o.048 



n /IQI +0.180+0.031 
U.^yi_Q^]_27-0.033 



0.855 



0.961 



+0.048+0.058 
-0.012-0.062 



-0.125+0.053 
-0.130-0.056 
1 f\f\o +0.051+0 mi 
l.UUO_o.074-0.065 



n fiSI +0.105+0.040 

U.UOi_Q QgQ_Q Q42 

n Qc;A+0-046+0.055 
"-'•^^^-0.022-0.059 
1 099+0-025+0.061 



order higher in logarithmic accuracy, (ii) instead of the scale choice /ig ~ ttih/N inherent in 
the moment-space formalism we set the scale as discussed in Section 4, and (iii) while the two 
formalisms are equivalent in the threshold region, the power-suppressed terms differ between 
the two formulations. Neither the additional higher-log contributions nor the scale-setting 
prescription can account for the difference. To compare the two scale-setting prescriptions, we 
have evaluated the effective-theory moment-space result (38) both with /ig ~ uih/N and with 
our choice of the soft scale and find that the difference is small. Using our default choice for the 
soft scale, the effective-theory moment-space result is 50.1 pb, very close to what is obtained 
in the traditional framework. The difference thus arises from power corrections suppressed by 
{1 — z) or 1/A^, respectively. If we use the same scale setting in the moment-space formula 
(38) and momentum-space expression (30), then the difference between the two formulations 
amounts to an overall factor of ^/z, see (39). To check that this factor indeed accounts for the 
difference, we have multiplied our momentum space formula (30) by ^/z. After adjusting the 
matching corrections, we find a = 49.9pb instead of 48.5pb.^ The factor ^/z appears artificial, 
since it does not occur in the fixed-order expressions. On the other hand, with this factor 
included, the singular terms are larger and amount numerically to 96% of the full NNLO 
result (without this factor, they amount to 86%). As stressed above, the threshold dominance 
is observed numerically but not enforced parametrically. For this reason, equivalent definitions 
of the leading contribution can lead to somewhat different results. 

To conclude our discussion, let us briefly discuss the case of Higgs production with a jet 
veto, i.e. the cross section for the production of the Higgs boson and QCD radiation with 
Pt = \pj^\ < Such a veto reduces the background to H W^W~ l^l^vv from 

tt-production with subsequent t M'^i/ decay. It was observed that the X-f actor for Higgs 

^Note that the enhancement is very counter-intuitive: the procedure amounts to multiplying the parton 
luminosity ffggir/z) by ^/z which is smaller than 1 over the entire integration range and nevertheless leads 
a larger result. The enhancement arises because the kernel is a distribution and its plus-distribution part is 
sensitive to the derivative of \/zffgg{T/z). 
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production gets reduced when such a cut is imposed [8, 58, 59]. For example, at the LHC 
with = 30 GeV and itlh = 165 GeV, reference [8] finds i^NLO ~ -^nnlo ~ 1-4, while the 
inclusive fT-factors are -RTnlo ~ 1-8 and i^NNLO ~ 2.2 (using MRST2004 PDFs). This implies 
that a significant portion of the perturbative corrections comes from the region pt > 30 GeV. 
Since the cut only excludes hard radiation, the threshold region should not be affected by it. 
At first sight, these large corrections from hard radiation seem difficult to reconcile with our 
finding that roughly 90% of the NLO and NNLO total cross sections arise from the partonic 
threshold terms. 

More formally, for a given value of z we have pt < ^^"2^^^ ' "^hich shows that the radiation 
from the threshold region z ^ 1 has small pt- The veto pt < p^^° therefore does not constrain 
the radiation above a certain value of 2; > zq. On the other hand, the converse is not true: 
small Pt does not imply that the radiation is soft. To isolate the amount of soft radiation 
present for a given pt cut, we have imposed a cut z > Zq in the integration over the leading 
singular terms. For pt = 30 GeV and Mh = 165 GeV, we have Zq = 0.7. We find that 
evaluating the threshold terms (7) with such a cut leaves them essentially unchanged. So we 
are forced to conclude that the cross section for pt < Pt*° is smaller than the contribution from 
the threshold terms alone, which implies that hard radiation with small pt, which comprises 
radiation with a small angle with respect to the beam, gives a sizable negative correction to 
the cross section. 

It is difficult to draw any conclusions from the above discussion, other than the trivial 
statement that also processes involving hard radiation can receive large higher-order correc- 
tions. To understand whether there is a physics reason for the observed reduction of the 
i^T-factors in the presence of a jet veto, it would be interesting to analyze Higgs production 
in association with a high-p^ jet in the effective theory. This process is mediated by opera- 
tors involving additional collinear fields for the partons inside the jet. The corresponding hard 
function will contain both space-like and time-like Sudakov double logarithms, which will need 
to be disentangled, before the enhanced contributions can be resummed and compared to the 
enhanced contributions affecting the total rate. 

6 RG-improvement for other time-like processes 

Having discussed Higgs-boson production in detail, we now brieffy explore the effect of choos- 
ing a time-like renormalization point /x^ < for other processes. Our treatment applies 
immediately to Drell-Yan production, but the numerical effects are less dramatic than for 
Higgs-boson production, as we pointed out in [18]. In addition, it is interesting to compare 
these production processes to inclusive decays such as e+e" hadrons, t ^ + hadrons, or 
hadronic Higgs-boson decay. For inclusive decay rates Sudakov double logarithms and the as- 
sociated TT^-enhanced terms are absent, since they cancel between real and virtual corrections 
by virtue of the KLN theorem [57, 60] . As a consequence, the effects of choosing //^ < are 
small unless the characteristic momentum scale is quite low. 
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Table 4: Predictions for the Drell-Yan cross section da /dq^ at \fs = 38.76 GeV for an invariant 
mass of ^/q^ — 8 GeV of the lepton pair. Units are pb/GeV^. 





fixed order 


threshold 


threshold + tt^ 


LO 
NLO 

NNLO 


U.Z»»_Q Q4Q 

n 440+0-051 

-0.041 

'-'■'-"-'■'-0.025 


n 436+0-062 
n 4qs +0-011 

n r,-| 9+0.002 
'-'■■'1- -0.001 


700+0-091 
p;p;q +0-014 

n r,o 1 +0.009 
1 _o.006 



6.1 Drell-Yan process 

Near the partonic threshold, the Drell-Yan cross section factors into a hard and a soft function, 
and threshold resummation proceeds in complete analogy to the Higgs case [28]. Instead of 
the scalar two-gluon operator (11), Drell-Yan production is mediated by the electromagnetic 
current q'~f^q. The hard function is given by the rcnormalized on-shell vector form factor 
Cv{—q^ 1 IJ^'ii) ■ The same hard function also appears in deep-inelastic scattering, but evaluated 
at space-like momentum transfer. The analytic continuation of the form factor to the time-like 
region produces tt^ terms, which were resummed in [19-22] . The Drell-Yan case can be treated 
in exactly the same way as Higgs-boson production. The quantity Cy fulfills a RG equation of 
the same form as (21) for C5, however the relevant cusp anomalous dimension F^^p is smaller 
by a factor Cp/Ca = 4/9. The resummation effects are thus smaller than in the case of Higgs- 
boson production and have the form exp(Ci?Q;s7r/2) at leading order, see (33). Comparing the 
expansion of the hard function at time-like and space-like renormalization points, we find 

|Cy(-9^9^)^ = 1 + 0.0845 + 0.0292 + ... , 

(44) 

|Cy(-9^ -q^)? = 1 - 0.1451 - 0.0012 + . . . . 

The two-loop correction is reduced for = — g^, however, at onc-Ioop order the correction 
increases since there is a partial cancellation between the vr^-enhanced terms and the constant 
piece for //^ > 0. Numerical results for the resummed Drell-Yan cross section are given in 
Table 4 for the case of E866/NuSea [61], i.e., proton-proton collisions at ^/s — 38.76 GeV 
and = (8GeV)^. The scales \fj,h\, yU-s, and /i/ in the resummation formula for Drell-Yan 
production have been chosen as in [28], and MRST2004 PDFs were used [52]. The numbers 
in the table include the matching to fixed-order perturbation theory at the corresponding 
order. At NNLO, the difference between ordinary threshold resummation with nl > and 
the combined resummation with //^ < is about 4%, significant because of the large ag value 
at such low energies. Convergence is similar in both cases, with negative instead of positive 
corrections for < 0. The numbers obtained in the two cases do not agree within their 
respective scale uncertainties. They would be compatible if the hard scale would be varied up 
and down by a factor 2 as we do in the present paper. However, a smaller variation < 
ijLh < 2>/g^ was used in [28]. In view of the disagreement, a variation by the conventional 
factor of 2 seems more appropriate. 

The main goal of the E866/NuSea experiment was to provide a determination of the sea- 
quark PDFs of the proton. The resummation of tt^ terms would affect this determination in the 
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form of an overall normalization factor. The absolute normalization of the cross section is also 
interesting at higher energies, e.g. for using the Drell-Yan process to monitor the luminosity 
at the LHC However, at NNLO the difference between the two scale-setting prescriptions 
scales as and would thus be four times smaller at = mz than at ^/q^ = 8 GeV, and 
completely negligible for higher-mass Drell-Yan pairs. 

6.2 e+e~ hadrons and hadronic r decays 

The total e^e^ hadrons cross section satisfies the relation 



where the sum extends over all quark flavors with 2mq < y/s, and the current-current vacuum 
correlator Ilqq{Q'^) is related to the Adler D-f unction as 



The quantity Ilqg{—q^ + ie) denotes its analytic continuation to the region of time-like mo- 
mentum transfer. For simplicity, we neglect the masses of the light quarks and assume that s 
is far away from quark thresholds. 

The Adler function in massless QCD is RG invariant, implying that its evolution equa- 
tion dD{Q'^) / din II = is trivially free of the cusp contributions associated with Sudakov 
logarithms and the large tt^ terms encountered in the case of Higgs-boson and Drell-Yan pro- 
duction. It follows that the perturbative expansion of the Adler function can be written as 

D(Q^) = l + ±dJ^)'', (47) 



n=l 



where the expansion coefficients dn are pure numbers, independent of the renormalization 
scale. Explicitly we have (setting Nc = 3) [62-64] 

di = 4, d2 = (22-16C3)/3o + ^~2.77/3o + 1.33. (48) 

In contrast to the case of Higgs-boson production, vr^-enhanced perturbative corrections enter 
the Adler function and the e^e~ hadrons cross section only at 0{al) and higher. Re- 
expressing the QCD coupling in terms of ^^(/x^), integrating relation (46), and analytically 
continuing to the time-like region, ^ —s + ie, we obtain 

Att' lmUqq{-s + ie)^l + d,^ + d, (^^^ ' + .... (49) 

This formula is routinely used in the calculation of the total cross section. Differences between 
the perturbative series in (47) and (49) arise starting at 0{al). 
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The RG-improved expression for the imaginary part of the correlator is obtained by ex- 
panding the result for the current-current correlator nyy(— s + ie) in powers of as{—s + ie), 
i.e., using perturbation theory at time- like momentum transfer [23, 65]. Integrating relation 
(46), we have 



4:7rllqq{—s + ie) ^ — 

TT 

= -ln(-s + ze) - — 

TT TT/io 



, D(a) 

da —rrT + const. 

/3(a) 

di In as{—s -\-ie) + I ^2 



diPi \ as{—s + ie) 



(50) 



/3o 



47r 



+ ... 



-I- const.. 



where f3i/ f3o = (102 — 38n//3)/(ll — 2n//3). Eliminating the time-like coupling in favor of 
the space-like one using (25), we obtain 



47r ImUqqi—s + ie) 



di 



as{s) arctan(as) 



47r 



+ 



An 



1 + 



diPi f arctan(as) ln(l -|- a^) 

^2/^0 \ 0.S 2 



+ ... , 

(51) 



where = (3oas{s)/A. This formula was first derived in [23]. It is the RG-improved version 
of (49), which should be used whenever is an 0{1) parameter. For example, with Ug — 0.7 
as appropriate for s — we find 



47r Imll 



qq\ 



-s + ie) = 1 + 0.872 d^ 



An 



+ d2 



aAs 



47r 



0.671 - 0.219 



d2Po 



+ ..., (52) 



which for nj = 4 yields a reduction of the two-loop coefficient by a factor of 0.45. For higher 
values of s corresponding to weak-scale processes, on the other hand, the modifications with 
respect to (49) are very small. 

Relation (51) shows that vr^-enhanced corrections to the e~^e~ — >■ hadrons cross sec- 
tions arise first at 0{al) in fixed-order perturbation theory. The leading term results from 
the expansion of the arctan(as)/a5 factor in the first line and yields (— 7r^/9Q/48)(as/7r)^ ^ 
— 14.3 {as/n)^ for Uf = 4, which is a rather large correction. For this reason, it was argued in 
[23] that for time-like processes as{s) = (4//9o) arctan(as) is a better expansion parameter than 
as{s). As an alternative choice, the authors of [65] suggested to use |q;s(— s)| ~ as{s)/ ^/T+~a^ 
for the expansion parameter. Note that both quantities have the property that they "freeze" 
in the infrared. 

Our RG-improved result (51) coincides with the expression for the cross section obtained 
using contour- improved perturbation theory [24, 25]. The analytic properties of the Adler 
function imply the relation 



\s'\=s 



d^D{e''^s). 



(53) 
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Inserting here the expansion (47) and using the generahzation 

of relation (25), where = {(p/tt) /3oQ;s(a*^)/4, we readily recover (51). The relevant contour 
integrals are evaluated in Appendix B. Weighted contour integrals over D functions can also 
be used to calculate the total hadronic decay rate of the r lepton and, more generally, moments 
of the corresponding spectral functions accessible in r decay [66, 67]. The equivalence of our 
approach with contour-improved perturbation theory extends to these cases as well. 



6.3 Hadronic Higgs-boson decay 

The total H ^ gg decay rate may be written as 

nH^gg) = ^^K. 



367r3V2 



2 ..2M2 TT 



(55) 



2m% 



lmUgg{-mH + ie, jj^) , 



where the correlator Tlgy is defined as in [68], and the leading-order 7^-factor is Ki^q = a^(yU^). 
We can use the RG-improved expression (20) for the matching coefficient Ct{m^, /x^) to rewrite 
the i^'-factor in the form 



K 



TT 



Ktimi) Knimij) , (56) 



where we have introduced the function 



bit,') = 



/?(a,(/.^))K(/.^) 



471 



+ ... . 



(57) 



The K-factor is independent of the scales /if and /i and factorizes into a product of two 
RG-invariant quantities depending on the scales rut and tuh- Evaluating the first factor with 
/If = TTit, we obtain (setting Nc = 3) 



X.(m?) = %^ . 1 + f 22 - ^ + 0{a',) . 



(58) 



Using the fact that the decay rate is RG invariant, i.e., that dK/dlnii = 0, the second factor 
may be written in analogy with (53) as 



TT 1 r ds 

Knimjj) = alipi') b'{ii') ^ Imn,,(-ml, + ie, /^') = ^ f ^ Dh{-s') . 



(59) 



\s'\=ml 



The perturbative series of the function Dh{Q') takes the form 



Dh{Q') = a^(g^) 



1 + ( 

n=l ^ 



(60) 
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where the expansion coefficients d„ are pure numbers, independent of the renormalization 
scale. At NLO we have [68] 



,H 70 14 

=73--n, + — 



(61) 



In fixed-order perturbation theory, the perturbative series for the function KH{m'\j) up to 
0{a^) takes the same form as (60) but with as{Q'^) replaced by Q;s(m|^). The corresponding 
expression in RG-improved perturbation theory is obtained by inserting the expansion (60) 
into (59) and using the contour integrals from Appendix B. At NLO we find 



KH{m\) 



al{my) 



1 + 



as{m]j) 



1 + { 1 + 47r 

where a = (3oas{mjj)/4:. For a = 0.2, we find numerically 
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An 



0.962 (if - 0.088^ 
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(62) 



(63) 



corresponding to a reduction of the LO and NLO terms by 4% and 8%, respectively. 

As a final comment, we emphasize again the different nature of the perturbative series for 
Higgs-boson production compared with that for Higgs-boson decay or the e~^e~ hadrons 
cross section. In the former case, vr^-enhanced perturbative corrections resulting from Sudakov 
double logarithms arise already at first order in perturbation theory and are due to the cusp 
anomalous dimension governing the RG evolution of the hard component of the hard-scattering 
kernel C. Terms proportional to n/, which reflect the evolution of the QCD coupling constant, 
enter first at two-loop order. The scale dependence of the ci;^(/i^) factor in the Born-level cross 
section is compensated by that of the gluon-gluon luminosity function. For the cases of Higgs- 
boson decay or the e+e" — > hadrons cross section, on the other hand, the perturbative series are 
free of 7r^-enhanced contributions up to NNLO, while n/-dependent corrections arise already at 
first non-trivial order in the expansion in ag- One may therefore argue that the resummation 
of terms of the form /^qCi"^^ is more important in these cases than the resummation of tt^- 
enhanced contributions. In fact, the leading such effects can be absorbed by setting the 
renormalization scale in the running coupling according to the BLM prescription [69]. The 
leading 0{Poal) correction in the series for the e+e~ — > hadrons cross section is absorbed 
by using the scale /xblm = ^'^(3-11/4, ^ ^ 0.708 ^/s in the running coupling constant. In the 
case of Higgs-boson decay one should use Hblm = c'^^^rriH ~ 0.174m//. The resummation of 
BLM terms to all orders and the associated renormalon ambiguities for the scalar correlator 
were analyzed in [70] , where also the resummation of the vr^-enhanced terms was discussed. 



7 Conclusions 

We have presented a RG-improved prediction for the total Higgs-boson production cross sec- 
tion. Our result is based on a factorization theorem for the partonic cross section near thresh- 
old, which is obtained by considering a sequence of effective theories in which the contributions 
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associated with higher scales are successively integrated out. In a first step the top-quark is 
removed and an effective Hgg interaction is derived. Subsequently, hard and soft modes are in- 
tegrated out and absorbed into Wilson coefficients H and S defined in SCET. This separation 
of the physics associated with different scales enables us to evaluate the different contributions 
at their optimal scale, where they have a well-behaved perturbative expansion. After this is 
done, the different elements are evolved to a common scale by solving the RG equations obeyed 
by the Wilson coefficients. 

The hard function H is given by the square of the on-shell gluon form factor, and we 
observe that this function receives large perturbative corrections due to tt^ terms arising in 
the analytic continuation of this form factor to time-like momentum transfer. These large 
corrections can be avoided by evaluating the hard function in terms of an expansion in powers 
of the running coupling at time-like momentum transfer. While the choice ^\ = —rri^j is 
unconventional, it is natural in the sense that it minimizes the logarithms in the hard function. 
The RG then sums up the logarithms that arise in the evolution between different scales. In 
our case, part of this evolution takes place in the complex /x^-plane and is driven by the cusp 
anomalous dimension, and the logarithms which are resummed are Sudakov double logarithms 
of —1. Our results for the resummed cross sections are stable under variation of the matching 
scales. Compared to the fixed-order expansion, the resummed perturbative series exhibits 
smaller theoretical uncertainties and better convergence. This improvement is mainly due to 
the resummation of the 7r^-enhanced terms, while soft-gluon resummation alone has a small 
impact. We find that the effects of RG improvement are significant even at NNLO. For 
example, for itlh = 120 GcV the resummed NNLO cross sections at the Tevatron and the 
LHC exceed the fixed-order predictions by about 13% and 8%, respectively. We perform a 
detailed analytical and numerical comparison with the traditional moment-space approach 
to resummation. While formally equivalent, the two methods differ by terms that are power 
suppressed in the threshold region. We analyze the difference in detail and show that it explains 
the somewhat larger threshold resummation effects found using the traditional method. 

We also discussed RG improvement for a number of other time-like processes, for which 
Sudakov logarithms are absent and our method reduces to contour-improved perturbation 
theory. The effects of choosing a time-like renormahzation scale are modest for the cases 
of the total e+e~ — > hadrons cross section and for the hadronic decay rate of the Higgs 
boson. Because these inclusive processes do not contain double-logarithmic corrections, the 
TT^-enhanced contributions arise solely from the running of and start at NNLO in fixed- 
order perturbation theory. These effects are significant for hadronic r decays, because the 
characteristic scale is rather low. For Higgs-boson decay, on the other hand, we have argued 
that BLM-type corrections of 0{l31^a'l+^) may be more important phenomenologically. 

It will be interesting to extend the methods developed here to other hard-scattering pro- 
cesses, such as jet production or heavy-quark production at hadron colliders. In general, 
processes containing collinear partons in both the initial and final states of the collision are 
characterized by an interplay of both time-like and space-like momentum transfers, so that 
the scale separation in the effective theory becomes more complicated. We hope to return to 
such issues in future work. 
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A Anomalous dimensions 

We write the perturbative expansions of the anomalous dimensions and /9-function in the form 

°° +1 °° +1 

7(«.) = E7n(g)" , PM--2a,J2^r^{^y . (Al) 

n=0 n=0 

The cusp anomalous dimension in the adjoint representation is given (at least up to three-loop 
order) by Ca/Cf times that in the fundamental representation. The first four expansion coef- 
ficients of the cusp anomalous dimension and the /3- function can be found, e.g., in Appendix B 
of [28] , where also the explicit expressions for the evolution functions S and op in (23) valid 
at NNLO are summarized. 

The first three expansion coefficients of the anomalous dimension 7*^ entering the evolution 
equation (21) of the hard matching coefficient Cs read [15, 71] 

7o =0, 

^ 160 llTT^ ^ \ / 208 An^\ 

7f = ^^1 ( + + 4(3) + C^TFUf (-—- — ]- SC^TFUf , 
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37045 61097r2 3197r^ /244 407r2 _ 



729 243 135 V 3 9 



16780 2 3 967r2 1647r^ 1424 . \ (^2) 



+ CiTpUf + + Cs 

^ V 729 243 135 27 ^ 

^ ^ ^ /1178 Att^ 16tt* 608 ^ \ ^^3^ 
+ CACFT^Uf i^— ^ Cs j + SC^TFUf 

^ _o 2 /24520 807r2 448 \ 176^ . o 

+ CaTIuI + Cs + CpTWf . 

■'^ V 729 81 27 ^7 9 ^ ^ 

The first three coefficients of the anomalous dimension 7^, which equals one half of the 
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coefficient of tfie 5{1 — x) term in tfie Altarelli-Parisi splitting function Pg^g{x), are [48] 



16 



7f = ( 3 + 3(3 1 - y CATprif - ACFTpUf , 
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79 47r2 llTT^ /536 Stt^ , ^ 



/233 Stt^ 27r4 160 ^ 
— CACpTpUf + 2CpTFnf + — CATpUf + — CpTpuj . 



(A3) 



B Contour integrals 

Here we evaluate the contour integrals 



1 r da' f a,(-s') Y 1 r /^ a,(e"'s) 

|S'|=S 



(Bl) 



where a = f3as{s)/4. Using relation (54) for the running coupling in the complex momentum 
plane, we find at NLO 



/n(a) = F{rj) 



77-»n-l f3o 47r 



+ 



r)—*n 



(B2) 



where 



Explicitly, we obtain for the first three integrals 



(B3) 



arctan(a) Pi as{s) 1 
Ii{aj = h 



/2(a) 



1 _^ /3i as(s) 1 



l + a2 /3o 47r (1 + a^)^ 
/3i q;5(s) 1 



+ 



+ Po 47r (l + a2) 



2^3 



^^^^-lln(l + a^)-l 

(1 - a^) ^^^^^ - ln(l + a^) - 1 
a 

(l-3a^)^^^^-^ln(l + a^)-l + 



a 
"3 



(B4) 
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